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Abstract 

It is shown that the dispersion and attenuation functions in a linear 
viscoelastic medium with a positive relaxation spectrum have a sublinear 
growth rate at very high frequencies. A local dispersion relation in para- 
metric form is found. The exact limit between attenuation growth rates 
compatible and incompatible with finite propagation speed is found. 

Incompatibility of superlinear frequency dependence of attenuation 
with finite speed of propagation and with the assumption of positive re- 
laxation spectrum is demonstrated. 

List of symbols. 

f{p) Jg°° f{t) e~P* dt Laplace transform; 

9{t) unit step function; 

t+ ^(01^1" homogeneous distribution. 
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1 Introduction. 



A very elegant theory linking attenuation and dispersion is presented for vis- 
coelastic media with positive relaxation spectrum. All the attenuation and 
dispersion functions compatible with the theory are represented by simple inte- 
gral representations. The attenuation function and the dispersion function are 
both expressed as transforms of a positive measure (the dispersion-attenuation 
measure). The measure is arbitrary except for a very mild growth condition. 
These expressions can be considered as a dispersion relation in parametric form. 

In contrast the acoustic Kramers-Kronig dispersion [2^ relations are non- 
local. They express the dispersion function in terms of the attenuation function 
or conversely. This presupposes that one of these functions (usually the disper- 
sion function) is very accurately known and consistent with the basic assump- 
tions of the theory. On the other hand, substituting any positive measure in the 
parametric dispersion relation yields an admissible (compatible with the theory) 
dispersion and attenuation. The origin of the Kramers-Kronig dispersion rela- 
tions is unclear. In electromagnetic theory they follow from the causality of the 
time-domain kernel representing the dielectric constant in the dispersive case. 
In acoustics they follow from an ad hoc assumption about the analytic prop- 
erties of the wave number. It is namely assumed that the wave number is the 
Fourier transform of a causal function or distribution. The physical meaning of 
the causal function or distribution is unclear hence the justification of the acous- 
tic Kramers-Kronig dispersion relation is missing. Various inequalities imposed 
on the complex continuation of the wave number cannot be expressed in terms 
of a constitutive assumption. Hence the acoustic Kramers-Kronig dispersion 
relation are an ad hoc addition to the constitutive equation, often incompatible 
with it. 

It will be also shown that in media with non-negative relaxation spectrum the 
frequency dependence of the attenuation function in the high frequency range is 
sublinear. In the case of power law attenuation the attenuation and dispersion 
are proportional to a power of frequency Numerous experiments in acous- 
tics indicate that the power-attenuation law accurately represents the frequency 
dependence of dispersion and attenuation over several decades of frequency. The 
theory based on positive relaxation spectrum implies that < a < 1. This seems 
to contradict the experimental ultrasound investigations of numerous materials 
which point to higher values of the exponent in the power law attenuation. For 
example in ultrasound investigations of soft tissues the exponent varies between 
1 and 1.5, while in some viscoelastic fluids such as castor oil it lies between 
1.5 and 2. We shall refer to this case as superlinear frequency dependence. 
Typical values of the power-law exponent in medical applications using ultra- 
sound transducers are a — 1.3 in bovine liver for 1-100 MHz, a = 1 — 2 in 
human myocardium and other bio-tissues [2^. a — 1.66 in castor oil at ca 
250 MHz. Values in the range 1-2 are observed at lower frequencies in aromatic 
polyurethanes |11) . Nearly linear frequency dependence of attenuation is well 
documented in seismology |10| . Approximately linear frequency dependence of 
attenuation has been observed in geological materials in the range 140 Hz to 
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2.2 MHz. 



Several papers have been devoted to a theoretical underpinning of the su- 
perlinear dispersion-attenuation models [25l [26l [6l |28l 13 El • In order to resolve 
some problems Chen and Holm [S] suggested to add a fractional Laplacian of 
order < y < 2 of the velocity field to the usual Laplacian of the displacement 
field in the equations of motion. Their paper still leads to an unbounded sound 
speed for y > 1 and adds a new problem: the equation of motion does not have 
the form of a viscoelastic equation of motion. 

Sublinear power law attenuation (a < 1) has also been reported in experi- 
mental investigations |21j . Non-power law attenuation laws are usually derived 
from the constitutive laws. Investigations of creep and relaxation in viscoelas- 
tic materials always support the assumption of positive relaxation spectrum 
(e.g. [El |9] for creep in metals, [17] for the upper mantle with a — 1/3) and 
therefore models derived from constitutive relations exhibit sublinear attenua- 
tion and dispersion at high frequencies. In particular this applies to the Cole- 
Cole, Havriliak-Negami and Cole-Davidson and Kohlrausch- Williams- Watts re- 
laxation laws commonly applied in phenomenological rock mechanics, polymer 
rheology [T], bio-tissue mechanics (e.g. for bone collagen [52]) as well as for 
ionic glasses [4]. 

Another abnormal feature of wave propagation in media with superlinear 
power laws (and more generally in media with superlinear asymptotic frequency 
dependence) is appearance of precursors. The precursors extend to infinity and 
thus the speed of propagation of disturbances is infinite. Finite speed of wave 
propagation requires that in the high-frequency range the exponent of the power 
law does not exceed 1. 

It is a very challenging problem how to explain the incompatibility between 
the theory and experiment in the superlinear case. It seems likely that the atten- 
uation observed at ultrasound frequencies significantly differs from the asymp- 
totic behavior of attenuation at the frequency tending to infinity. 

One might surmise that the frequency range of ultrasound measurements is 
still far below the asymptotic high frequency range in which different mecha- 
nisms are at play. This suggests studying models of attenuation with a slowly 
varying power-law exponent. 

2 Constitutive assumptions and basic definitions. 

In viscoelasticity the relaxation modulus G, defined by the constitutive stress- 
strain relation 



Jo 

is assumed to have positive relaxation spectral density h. The latter statement 
means that for every i > 




(1) 




t > 



(2) 
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where r — l/r is the inverse of the relaxation time and h{r) > 0. Eq. (|2| 
represents a superposition of a continuum of Debye elements. For mathematical 
convenience eq. ([2| will be replaced by a more general equation 

G{t) = [ e-*Xdr), t>0 (3) 

J[0,oo [ 

where ^ is a positive measure: fi{[a,b]) > for every interval [a,b] of the pos- 
itive real axis. As indicated in the subscript of the integral sign the range of 
integration is the set of reals satisfying the inequality < r < c». In general the 
measure /i({0}) of the one-point set {0} is finite and equal to the equilibrium 
modulus Goo ■— limt_i.oo G(t). An additional assumption 

— - < oo (4) 

[o,oo[ ^ + r 

ensures that G is integrable over [0, 1], [13]. The relaxation modulus assumes a 
finite value G(0) = M at if the measure fi has a finite mass M. 

The right-hand side of eq. ([3) can be replaced by a Stieltjes integral with 
respect to the function g{r) = fi( 0,r]): 

G{t)^ f e-'^dg{r) (5) 

J [0,00 [ 

The function g is non-decreasing and right-continuous: g{r) = /i([0,r]) = 
lime_).o+ /^([O, r + e]), and g{r) = for r < 0. 

In contrast to eq. ^ the integral representations ^ and ([s]) include as 
special cases finite spectra of relaxation times corresponding to superpositions 
of a finite number of Debye elements 

N 

G(i) = ^c„e-'^"*, c„>0,r„>0 forn = l,...,iV (6) 

n=l 

(Prony sums), infinite discrete spectra (corresponding to Dirichlet series) as 
well as discrete spectra embedded in continuous spectra. However the main 
advantage of ([3| over ([2| is the availability of very powerful mathematical theory 
which ensures logical equivalence of certain statements about material response 
functions. 

In particular a function satisfying (|3| is completely monotone. A function G 
is said to be completely monotone if it continuously differentiable to any order 
and 

(-1)" D" G{t) > for all non-negative integers n and t > (7) 

Bernstein's theorem asserts that eq. ([s]) is equivalent to J7|. There is no 

such simple characterization of functions which have the form H. 

It was established in [3l [181 US] that viscoelastic relaxation moduli are com- 
pletely monotone. 
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For us the main benefit from using ([3| instead of ([2| is the equivalence of (|3| 
with a property of the dispersion and attenuation that will be explained below. 
That is, certain statements about attenuation and dispersion functions follow 
from Q or, conversely, imply that ([S]) does not hold. 

We now recall the definition of a Bernstein function 2J. A differentiable 
function / on ]0, (X)[ is a Bernstein function if / > and its derivative /' 
is completely monotone. A Bernstein function is non-negative, continuous on 
]0, oo[ and non-decreasing, hence it has a finite value at 0. A function / on 
[0, oo[ which has the form f{x) — g{x) for some Bernstein function g is called 
a complete Bernstein function (CBF) |15j . It can be proved that every complete 
Bernstein function is a Bernstein function [15| . 

The following facts about complete Bernstein functions will be needed here. 

Theorem 2.1 

A real function f on ]0, cx)[ is a complete Bernstein function if (i) f has an 
analytic continuation f{z) to the complex plane cut along the negative real 
axis; (ii) /(O) > 0, (Hi) f (z) = ,f{z), (iv) Im/(z) > in the upper half plane 
Imz > 0. 



Theorem 2.1 has the following corollaries: (1) If / is a CBF then /" is a CBF 
if < a < 1 d4J. 

(2) If / ^ then the function x/f{x) is a CBF. 

Every complete Bernstein function / has the integral representation: 

f{x) = a + bx + x f (8) 
■/]o,oo[ x + r 

where a,b > and v is a positive measure satisfying the inequality 

TT7 < (9) 

[0,oo[ 

m- 

If G is a completely monotone function integrable over [0, 1] and satisfying 
eq. ((sl) then 

G(p) = / ^ = + [ ^ 

It follows that 

Qb):=pG(p)=M{0})+ / ^ (10) 
is a complete Bernstein function [l4] . 



3 Dispersion and attenuation. 

The Green's function Q{t, x) is defined as the solution of the problem 

puM ^ G{t) * u^txx + S{x) S{t) (11) 
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with zero initial data. In a three-dimensional space 



where 



-1 



(27T)3 



ioo+e 



Q{p) 



dp 



ioo-fe 



fcdfc 



87t2i 



— loo+e 



Q{p) 



pPt-B(p)r 



B{p) -.^p 



dp (12) 



(13) 



Q(p)i/2 

If G is completely monotone then Q is a complete Bernstein function. If a 
function / is a complete Bernstein function then the functions /^/^ and p/ f{p) 
are complete Bernstein functions. Hence _B is a complete Bernstein function 
and B{p) = a + cp + b{p) where a, c > and the dispersion-attenuation function 
b has the integral representation 



Hp) -^p 

J]0,c 



v{dr) 



(14) 



and 1/ is a positive measure satisfying A function b having an integral 

representation (14) will be called an admissible dispersion- attenuation function. 



The measure v is the spectral measure of the dispersion-attenuation function 
b. A dispersion-attenuation function 6 of a viscoelastic material with a positive 
relaxation spectrum is admissible. Conversely, for any admissible dispersion- 
attenuation function there is a completely monotone relaxation modulus G sat- 



isfying eq. ( 13 I and eq. ( 10 1 



Using this equivalence it is possible to decouple the study of admissible 
dispersion-attenuation functions from considering specific constitutive equa- 
tions. Furthermore, in order to avoid incompatibility with viscoelastic theory 
experimental studies should target the spectral measure v oi b rather than the 
dispersion-attenuation function b{p). Any positive measure v satisfying ([9| is 
compatible with the theory and the implications of the theory for the dispersion- 



attenuation function b are encapsulated in (14). Using the spectral measure v 



instead of b to match the experimental data is analogous to expressing relax- 
ation modulus in terms of the relaxation spectrum, e.g. by applying Prony sums 
to represent experimental data. A discretization of the dispersion-attenuation 
measure is presented in Sec. |5.1| 

Let Goo = hmi_j.oo G{t). Since limp^o pG{p) — Goo, the inequality Goo > 
and eq. ( 13 ) imply that B{0) = and therefore a = 0. The limit limp_j.oo B{p) jp 



c implies that c = I/cq, where cq denotes the wave front speed, 
elusion follows from the fact that for p > 1 



The last con- 



1 

p + r 



< 



1 
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and the right-hand side is integrable with respect to the measure v in view of 
eq. ([9]) . For p — > cxi the integrand of 



]o,oo[ P + r 

tends to zero, hence, by the Lebesgue Dominated Convergence Theorem, the 
integral tends to zero as well. It follows additionally that b{p) ~ o[p] for p — > 00. 
Hence ( 12 1 can be recast in a more explicit form 

1 Aca+e 1 

g{t, x) - / — eP (*-I-I/-o)-Kp) r dp (15) 

and the dispersion-attenuation function h{p) has sublinear growth. 
We also note that the attenuation function is non-negative: 

A(j>) := Re h[p) = [ bl' +^ ^^P u{dr) > (16) 

7)0,00 [ \P + 1T 

for Rep > 0. 

The derivative 

b'ip) = b{p)/p - I < b{p)lp (17) 

Moreover 

i]o,oo[ (P + ~ 

Hence b'{p) = o[l] for p — > 00. 

The attenuation function A and the dispersion function I?(|)) := — Im6(p) 
satisfy linear dispersion equations in parametric form: 

A/ \ f bp + ^ R-ep , , , 



P(p) = - Imp I 

J]0,c 



r 



\p + r\ 



-,v{dr) (19) 



The measure v represents the dispersion-attenuation spectrum. An elemen- 
tary dispersion-attenuation is represented by the function p/ [p + r) for a fixed 
value of r. 



4 Implications of finite speed of wave propaga- 
tion. 

Since \/Q{p) = pJ{p) — J{0) + J'{p), where J is the creep compliance, we can 
express the Green's function ( |15[ ) in the form of a convolution 

g{t,x) = J'it) * H{t - \x\/cQ, \x\) + Ja Hit - |x|/co, |x|) 
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where H{t, r) is the inverse Laplace transform of the function e ^^^'^ (inr). In 
terms of the inverse Fourier transformation 



HiT,r) 



1 



dw 



Note that 6(— iw) — b{iuj). The integrand is square integrable if r > and 



(20) 



If eq. (20) holds then the Paley- Wiener theorem (Theorem XII in [TH]) can 
be applied. By this theorem H{t, r) = for r < if and only if 



Re6(-ia;) 
1+^2 



duj < oo 



(21) 



J' is a causal function. Hence, if H(t, r) = for t < then 

ft-\x\/co 



Q{t,x) 



vanishes for t < \x\/cq. 



J'(s)H{t- |a;|/co-s,|a;|)ds 



In particular, if b{p) ap", a > 0, Rep > 0, eqs (20) and (21 1 are ensured 
by the inequality < a < 1. A precursor appears for a > 1, as can be seen 
from Fig. [l] For a = 3/^ there is no wavefront and the peak is preceded by 
a precursor extending to infinity. The limit case, as we already know, is the 



^For a = 1/3, 2/3 the ct-stable probability can be expressed in terms of Airy functions, see 



mi- 
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asymptotic behavior b{p) and it entails unbounded propagation 

speed. 

For a general viscoelastic medium with a positive relaxation spectrum we 
note that 

ReK-iu.)=c.Im/ ^ = / (22) 

If the total mass M of v is finite then 

lim Re6(-iu;) = M (23) 

by the Lebesgue Dominated Convergence Theorem. In the general case the 
inequality 

CO 2uj 2uj 

valid for co > I and the Lebesgue Dominated Convergence Theorem imply that 
Re6(— iw) — o[uj] at w — > oo. 

We have thus proved that the asymptotic growth of Re b{—iuj) in the high- 
frequency range is sublinear. This is however insufficient to ensure convergence 



of the left-hand side of (21 1 and vanishing of the wave field ahead of the wave- 



front = Cot. For example for 6(— iw) ~oo w/ln"(a;) the integral in (21) does 
not always converge. Since we are interested in the high-frequency behavior, we 
can replace 1 -t- by cj^ in the denominator of (21 1. For a ~ 1 



y My) 



is unbounded for uj — > oo, hence (21) is not satisfied. For a ^ 1 

dy ln^~"(aj) — 1 

y ln"(y) " 

is unbounded for a < 1 and bounded for q > 1. We thus see that in contrast to 
the power law attenuation, linear growth is not a limit case for finite propagation 
speed. Some sublinear cases will also exhibit precursors ahead of the wave front. 

In particular the function b{p) = p/ln"(l + p) with < a < 1 is a CBF 
with the asymptotic properties discussed in the previous paragraph. Indeed, 
for Imp > the argument V' = arg(l -|- p) satisfies the inequality < ^ < n. 
Hence ln(l + p) maps the upper half plane into itself and is non-negative for 
p > 0. Hence, by Theorem |2.1[ ln(l -I- p) is a CBF. The same is true for 
ln"(l -l-p) if < a < 1. In view of a property of CBF functions mentioned after 



Theorem 2.1 this implies that b{p) is a CBF if < a < 1. Moreover 6(0) = and 
limp_>.oo b{p)/p = 0, hence b{p) is an admissible dispersion-attenuation function. 
We thus have produced an example of an admissible dispersion-attenuation 
function b{p) such that 6(— iw) has sublinear growth in the high frequency range 



but (21) is not satisfied. 
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5 Examples. 



5.1 Discrete dispersion-attenuation spectra. 

If the measure v = Cn£r„, fn-,Cn > 0, where Eq = 6{r — a) denotes the 

Dirac measure concentrated at the point a, then 

N 

=p^c„/(p + r„) (24) 

n=l 



Eq. ( 24 ) can be used to construct numerical approximations of experimental 



data. Using the relation Q{p) B{p)^ — pp^ and the assumption p ~ X]m=i '^m 
with dm > 0, Sm > for 1 < m < M , Sm ^ fn for all n, m, it is possible to 
express the relaxation data d„nSm hi terms of the dispersion-attenuation data 
c„ , r„ and conversely. 

5.2 Power-law attenuation. 

Power-law attenuation is commonly used to match experimental dispersion and 
attenuation data for a wide variety of real viscoelastic materials such as poly- 
mers, bio-tissues and some viscoelastic fluids. 

Consider the viscoelastic medium defined by the following equation of motion 

p (D^ + 2a Dj,+" + D^") u = A u + 5{t) 5{x) (25) 

with the initial data m(0,x) = u^t(0, a;) = in d dimensions, d = 1,3. Dg 
denotes the Caputo fractional derivative of order a 

Dg/(x):= f^l-^fis)ds (26) 
Jo r(l-a) 



It is assumed that ^ > 0, a > 0, < a < 1. The Laplace transformation 
Ct with respect to the time variable and the Fourier transformation J-^ with 



respect to the spatial variable bring eq. (251 to the following form 

pg{pfu{p,k)^-Ak^u{p,-k) + l (27) 



where u := Tx {Ct{u)) and g{p) := p -\- ap". Eq. ([27| implies that in the power 
law attenuation model Q{p) = p A/g{p)'^ and B{p) = g(p)/cQ. 

We shall begin with solving eq. ( [25| in one dimension. Applying the inverse 
Fourier transformation to eq. (|27|): 



1 f°° Qikx+pt 

27Ty„oo pgipy + Ak^ 

The contour can be closed by a large half-circle in upper-half complex fc-plane 
if X > 0, in the lower- half complex fc-plane if x < and 



1 1 



Ak^ + pg{p)^ 2ig{p)A 



1 1 



k-ig{p)lcQ k + ig{p)/co 
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where cq '■— y^A/p. We now restrict ourselves to imaginary values oi p — ~iw. 
Since Im[i (;(— iw)] = Keg{~iw) = a|w|" cos((l — a)n/2) > 0, the residuum at 
±ig{p)/co contributes if ±x > 0. Hence 

I Q-a{p)M/co+pt 

= ^ 9ip) 

and 

U^^\t,x) = -— QPi-9{p)M/co^ (28) 

AmAj^,^g{p) 



The solution u'^) of ^ in a three-dimensional space is given by the formula 



Hence 



1 /■ 

87T^i|a;|coA ' 



47tai/" Co A 

where 



1C30 

^ P„((t-N/co)/Hx|)") (29) 



27ti 



lOO 



Po.{z):^^i ey^e-y^'dy (30) 



lOO 



The function Pa is a totally skewed Levy stable probability density [231 E2] ■ 

In the case of power-law attenuation Re6(— iw) = aw" cos(7Ta/2) > and 
for p = —iLj 



QPt-g{p) ■r/co 



^— aw" cos(Q7r/2) 



Hence ([3U], Sec. 4.44) the integrals 



fit.r) / (-!)>-" p"g(p)"+ief '■/^Mp 



are uniformly convergent for r > and the derivatives S""*""* c^r™ exist 
for all positive integers n, m and r > 0. 

All the properties of the dispersion-attenuation function b were derived from 
the fact that B{p) — pp^^^/ Q(p)^/^ is a complete Bernstein function. It follows 
from the above property of B that Q(p)^/^ is a complete Bernstein function 
[Mj. It does not follow from here that Q is a also complete Bernstein function. 
Therefore G need not be a completely monotone function. A counterexample 
is provided by the power law attenuation model. If < a < 1/2 then the 
relaxation modulus G is not completely monotone but the attenuation function 



satisfies eq. (14 1. The latter condition is satisfied for the power law attenuation 



if and only if < a< 1/2 (Fig. m 
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Figure 2: The relaxation modulus G for the power-law attenuation. The expo- 
nent values are a — 0.3, 0.4, 0.5, 0.6, 0.9, from bottom to top. 

Theorem 5.1 

Tiie spectral measure of the dispersion-attenuation function ap°' , < a < 1, is 



7T 



Proof. The Laplace transforms 

/"OO 

a;"-i = / e~-y [y-"/r(l - a)] dy 
Jo 

y-" = / e-^^ [z"-Vr(a)] dz 
Jo 



hence 



_ sin(Q;7t) / ~ z" 



X + X 



dz 



which implies eq. (31 1 



(31) 

(32) 
(33) 



□ 



6 Superlinear power law attenuation. 

A large number of papers have been devoted to the implications of the Kramers- 
Kronig dispersion relations for the wave number K{lu). The Kramers-Kronig 
dispersion relations would follow from the assumption that the function K{u;) — 
uj /co is the Fourier transform of a causal function or causal tempered distribution 
L{t). A priori the function L{t) has no physical meaning and the assumption 
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of causality of L is unwarranted. Causality of L would however be justified by 
the assumption that 

B{p)=plco + L[p) 



and the equation of motion has the following form 

Cq"^ u,tt + (2/co) L * u,t + L * L * u 



(34) 



In this case iir(w) = i_B(— iw) = oj / co + \L{—\uj) =iw/co + i(w). Eq. (34) is how- 
ever incompatible with the viscoelastic constitutive equation. In a viscoelastic 
equation of motion integral operators should act on the Laplacian of u. 

In the authors try to guess the viscoelastic constitutive equation compat- 
ible with ( |34[ ) . Their approach involves an approximation of a spatial derivative 
by a temporal derivative. It is however possible to avoid an approximation by 



shifting the dispersive terms on the left-hand side of ( 34 ) to the right-hand side 



The Laplace transform of the left-hand side of eq. ( 34 1 is 



l + CoL{p)/p /cj 



assuming that u{0,x) — u^t{0,x) — 0. Hence (34) has the form 

Co~^u.tt = [G(t) * u^tx]^^ 
where G is the inverse Laplace transform of 

1 + Co L{p)/p 



(35) 



Expression ( 35 ) is the Laplace transform of a completely monotone function if 

r — ~ 1 ^ 

and only if 1 -I- cq L{p) jp /p \s the Laplace transform of a Bernstein function 

/ mm- 

For L{p) — ap"^, with a, a > 0, the function / assumes the form 6{t) + 
2coat+""/r(l-a)-hc§a2t+^^""Vr(3-2a). It is obvious that / is a Bernstein 
function (and G is completely monotone) if and only if 1/2 < a < 1. We also 



note that L(t) 



at. 



/r(— a) is a distribution. Convolution with L{t) is 



a Riemann-Liouville fractional differential operator of order a. The order of 



the fractional differential equation (34) is 2 if a < 1 and 2a if a > 1. (In the 



literature the highest-order derivative L * L * u is often incorrectly neglected). 
For a > 1 the order of the time-differential operator exceeds the order of the 
spatial differential operator, hence the equation is formally parabolic. 

Concerning superlinear power law attenuation, note that the logarithmic 
attenuation rate 



A{lo) := ReB{-iuj) = a cos(a7t/2) 



(36) 



is assumed non-negative, hence 1 < a < 3 must entail that a < 0. As the 
frequency tends to infinity the frequency-dependent phase speed c(a;), given by 
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the equation l/c(w) := Re[i-B(— iw)] = l/cg +a|a;|"~-^ sin(Q!7T/2) increases to 
its maximum value cq if < a < 1 ("abnormal dispersion"). For 1 < a < 2 
it increases from cq at zero frequency to infinity at a finite frequency uJi = 
[cq \a\ sin(a7T/2)]^^^"~"^'' and changes sign. This behavior is clearly unphysical. 
For 2 < a < 3 phase speed decreases from cg at cj = to at infinite frequency 
("normal dispersion"). 



7 Models with variable attenuation exponent. 

The limits on the dissipation-attenuation exponent a actually apply to the 
asymptotic value of b{p) at infinity. Is it possible that the experimentally mea- 
sured power law behavior applies to the middle frequency range? 

Define the variable dissipation-attenuation exponent as the function 

a(p)=ln(&(p))/ln(p), p>l (37) 

so that b{p) = p"^P\ This definition has a major flaw: a singularity at p = 1. 
The exponent decreases to — oo at p — >■ 1— and restarts from oo to decrease 
towards its asymptotic values. 

The simplest examples of dispersion-attenuation functions with variable ex- 
ponent are 

6i(p) =cK/(a+p"), 0<a<l (38) 

and 

h2{p)^c{l + Tp)°'-'^ [tp)^, 0</3<a<l (39) 

hi is a complete Bernstein function because p/bi{p) = p/c+ap^^°' /c is obviously 
a complete Bernstein function. Moreover foi(O) = and limp_j.oo &i (p) /p = 0. 
Hence bi is an admissible dispersion-attenuation function. 

In ord er to prove that 62 is a complete Bernstein function we need Thc- 
We now note that arg62(p) = arg(l -I- p)"~^ + argp^ < (a — 



2.1 



orem 

/3) a,Tgp + /3 argp — a argp. For p in the upper half-plane < argp < tt, hence 
< arg62(p) < 7T and 62 is a complete Bernstein function. Since 62(0) = and 
limp_>.tx) bi{p)/p = 0, 62 is admissible as a dispersion-relaxation function. In the 
first case, b — bi, the attenuation exponent a{p) < and a(p) from below 
as p — 00. In the second case 6 = 62 and a{p) increases from a(0) = /3 to 
a{oo) — a. It is thus likely that the exponent assumes a larger value in the high 
frequency range. Thus a value of the exponent above 1 in the middle frequency 
range is not likely. 



8 Conclusions. 

The class of admissible dispersion and attenuation functions can be charac- 
terized by a class of Radon measures. The theory applies only to sublinear 
attenuation growth in the high frequency range. 
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Supcrlincar growth of the attenuation function in the high frequency range 
is incompatible with the assumption of positive relaxation spectrum underlying 
theoretical and experimental viscoelasticity. Superlinear growth of attenuation 
also implies that the phase speed is unbounded for high frequencies and the 
main signal is preceded by a precursor of infinite extent. 

Explanation of superlinear frequency dependence of attenuation observed in 
many real materials remains a challenging task. 

References 

[1] R. L. Bagley and P. J. Torvik. On the fractional calculus model of vis- 
coelastic behavior. J. of Rheology, 30:133-155, 1986. 

[2] C. Berg and G. Forst. Potential Theory on Locally Compact Abelian Groups. 
Springer- Verlag, Berlin, 1975. 

[3] D. R. Bland. The Theory of Linear Viscoelasticity. Pergamon Press, Ox- 
ford, 1960. 

[4] G. Carini, M. Cutroni, M. Federico, G. Galli, and G. Tripodo. Phys. Rev. 
B, 30:7219-, 1984. 

[5] W. Chen and S. Holm. Fractional Laplacian time-space models for linear 
and nonlinear lossy media exhibiting arbitrary frequency power-law depen- 
dency. J. Acoust. Soc. Am., 114:2570-2574, 2003. 

[6] W. Chen and S. Holm. Modified Szabo's wave equation models for lossy 
media obeying frequency power law. J. Acoust. Soc. Am., 114:2570-2574, 
2003. 

[7] R. S. C. Cobbold, N. V. Sushilov, and A. C. Wcathcrmon. Transient prop- 
agation in media with classical or power-law loss. J. Acoust. Soc. Am., 
116:3294-3302, 2004. 

[8] E. N. da Andrade. On the viscous flow of metals and allied phenomena. 
Proc. Roy. Soc. London, A84:l-12, 1910. 

[9] E. N. da Andrade. On the validity of the t^^^ law of flow of metals. Phil. 
Mag., 7:(84), 1912. 

[10] W. Futterman. Dispersive body waves. J. Geophys. Res., 67:5279-5291, 
1962. 

[11] J. F. Guess and J. S. Campbell. Acoustic properties of some biocompatible 
polymers at body temperature. Ultrasound Med. Biol., 21:273-277, 1995. 

[12] A. Hanyga and M. Seredyhska. Asymptotic wavefront expansions in heredi- 
tary media with singular memory kernels. Quart. Appl. Math., LX:213-244, 
2002. 



15 



[13] A. Hanyga and M. Seredynska. Relations between relaxation modulus and 
creep compliance in anisotropic linear viscoelasticity. J. of Elasticity, 88:41- 
61, 2007. 

[14] A. Hanyga and M. Seredynska. Positivity of viscoelastic Green's functions. 
\arXiv:0908.307^ 1 [cond-mat], 2009. 

[15] N. Jacob. Pseudo-Differential Operators and Markov Processes, volume I. 
Imperial College Press, London, 2001. 

[16] J. F. Kelly, R. J. McGough, and M. M. Meerschaert. Analytical time- 
domain Green's functions for power-law media. J. Acoust. Soc. Am., 
124:2861-2872, 2008. 

[17] J. B. Minster and D. L. Anderson. A model of dislocation-controlled rhe- 
ology for the mantle. Phil. Trans. R. Soc. London, 299:319-356, 1981. 
Andrade creep. 

[18] A. Molinari. Viscoelasticite lineaire and functions completement mono- 
tones. J. de mecanique, 12:541-553, 1975. 

[19] R. E. A. C. Paley and N. Wiener. Fourier Transforms in the Complex 
Domain. AMS, New York, 1934. 

[20] I. Podlubny. Fractional Differential Equations. Academic Press, San Diego, 
1998. 

[21] A. Ribodetti and A. Hanyga. Some effects of the memory kernel singular- 
ity on wave propagation and inversion in poroelastic media, II: Inversion. 
Geophys. J. Int., 158:426-442, 2004. 

[22] N. Sasaki, Y. Nakayama, M. Yoshikawa, and A. Enyo. Stress relaxation 
function of bone and bone collagen. J. of Biomechanics, 26:1369-1376, 
1993. 

[23] Ken-Iti Sato. Levy processes and Infinitely Divisible Distributions. Cam- 
bridge University Press, Cambridge, 1999. 

[24] M. Seredynska and A. Hanyga. Cones of material response functions in ID 
and anisotropic linear viscoelasticity. Proc. Roy. Soc. London A, 465:3751— 
3770, 2009. 

[25] T. L. Szabo. Time domain wave equations for lossy media obeying a fre- 
quency power law. J. Acoust. Soc. Am., 96:491-500, 1994. 

[26] T. L. Szabo. Causal theories and data for acoustic attenuation obeying a 
frequency power law. J. Acoust. Soc. Am., 97:14-24, 1995. 

[27] T. L. Szabo and J. Wu. A model for longitudinal and shear wave propaga- 
tion in viscoelastic media. /. Acoust. Soc. Am., 107:2437-2446, 2000. 



16 



[28] K. R. Waters, M. S. Hughes, G. H. Brandenburger, and J. G. Miller. On a 
time-domain representation of the KramersKronig dispersion relations. J. 
Acoust. Soc. Am., 108:2114-2119, 2000. 



[29] R. L. Weaver and Y.-H. Pao. Dispersion relations for linear wave propa- 
gation in homogeneous and inhomogeneous media. J. Acoust. Soc. Am., 

22:1909-1918, 1981. 

[30] E. T. Whittaker and G. N. Watson. A Course of Modem Analysis. Cam- 
bridge University Press, Cambridge, 1927. 4-th ed. 

[31] D. V. Widder. The Laplace Transform. Princeton University Press, Prince- 
ton, 1946. 

[32] V. M. Zolotarev. Modem Theory of Summation of Random Variables. VSP, 
Utrecht, 1997. 



17 



